Mon. Not. R. Astron. Soc. 000,[T]l5](2008) Printed 16 September 2008 (MN I^TbX style file v2.2) 



Puzzling X-rays from the new colliding wind binary WR 65 (WC9d) 
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ABSTRACT 

We report the discovery of variability in the X-ray emission from the Wolf-Rayet type star 
WR65. Using archival Chandra data spanning over Syr we detect changes of the X-ray 
flux by a factor of 3 accompanied by changes in the X-ray spectra. We believe that this X- 
ray emission originates from wind-wind collision in a massive binary system. The observed 
changes can be explained by the variations in the emission measure of the hot plasma, and by 
the different absorption column along the binary orbit. The X-ray spectra of WR 65 display 
prominent emission features at wavelengths corresponding to the lines of strongly ionized 
Fe, Ca, Ar, S, Si, and Mg. WR65 is a carbon rich WC9d star that is a persistent dust maker. 
This is the first investigation of any X-ray spectrum for a star of this spectral type. There are 
indications that the dust and the complex geometry of the colliding wind region are pivotal in 
explaining the X-ray properties of WR 65. 
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1 INTRODUCTION 

More than 54% of all Wolf-Rayet (WR) star s are in binary sys- 
tems that consist of a WR and an OB-type star ( IWallacel2007n . The 
collision between two stellar winds in these systems produces char- 
acteristic signatures in differe nt wavelength bands th at can include 
non-thermal radio emission (lEichler&Usovlll993h , copious and 
variable X-rays l lStevens. Blondin & Pollock" 1992"), and dust emis- 
sion in the infrared (IR) iwilliams. 7a.n der Hucht & Thelll987l) . 

Subject of this Letter is the WC9d star WR65. 

Iwilliams. van der Hucht & Thd ( Il987h reported the presence 
of warm dust in this object. They discussed that in early-type 
WC stars an increase in the wind density provided by e.g. shock 
compression in colliding wi nds in binary systems may be sufficient 
to result in grain formation. IZubkol ( 1 1 998h found no need to invoke 
binarity to explain the presence of dust in late WC stars. He fitted 
the observed IR spectrum of WR 65 assuming a single star with 
stellar mass-loss in the form of dust with = 5.4 x 10"'" M© yr"' 
and surrounded by a dust shell with a n inner radius of 520^.. 
iLeitherer. Chapman & Koribalskil l ll997h detected radio emission 
from WR 65, but could not constrain its spectral index. 

In a previous paper ( Oskinova et al. 2003j) we demonstrated 
that single WC type stars are X-ray quiet and proposed that all X- 
ray active WC stars must be in binary systems. X-ray emission from 
WR 65 was detected and tentatively explained as a result of wind- 
wind collision. In this paper we report on the X-ray light curve 
of WR 65 and changes in its X-ray spectra that unambiguously 
confirm its status as a colliding-wind binary (CWB). 
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Table 1. Coordinates of WR65 from different sources 



Ref 


Band 


RA J2000 


DEC J2000 


Simbad 


opt. 


ISh 


13™ 41^:68 


-59° 11'43'.'3 


Chaoman et al. (19991 


radio 


15" 


13"- 41^49 


-59° 11'43'.'8 


Chandra HRC-I 


X-ray 


15'' 


13""41':76 


-59° ll'44'.'l 


Catalog USNO-Bl.O 


opt. 


15'- 


13"" 41^73 


-59° 11'43'.'4 




Figure 1. Archival Spitzer image (8 pim). WR 65 is the bright object at the 
bottom, close to the detector edge. The coordinates are galactic. 



2 WR 65 AND ITS X-RAY LIGHT-CURVE 

The coordinates of WR 65 from different sources are compiled in 
Table[T] The coordinates from Simbad and the USNO-Bl.O cata- 
log agree well with the source location in the Chandra HRI-I im- 
age. Except for WR65 there are no other known objects at this 
position. WR65 is a probable member of the cluster Pismis20, 
which contains one more WR star, WR67, and has a distance of 
d = 3272 ± 303 pc (Q'urner 1996). 

We fit the optical spectrum of WR 65 and its available pho- 
tometry using the Potsdam Wolf-Ra yet (PoWR) stellar atmosphere 
models jHamann & Grafeneill2003l) . The deduced stellar parame- 
ters are listed in Table|2] 
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Table 2. Stellar parameters of WR65 (from Hamann et al., in prep.) 





logLbol/io" 


ii„ [kms '] 


M[Moyr"']^ 


2.4 


5. 8. ..6.4 


» 2000 


lO"* " ... lO"* ' 



(a) depending on the adopted reddening law 

(b) assuming a clumping contrast Z) = 10 




Figure 2. Archival Chandra ACIS-I image (0.2-1 2.0 keV) of WR65. The 
observation is from 18 Oct. 2005. The bright object in the upper-left corner 
of the image is the pulsai' PSR B 1509-58. WR 65 is projected onto the out- 
skirts of the pulsai' wind emission. The coordinates are equatorial (J2000). 
North up, east left. 



A wind-blown bubble with a radius of 7' i s known around 
WR 65 ( lMarstonlll997h . lGiacani & Dubneil | |2004|) discovered a H i 
shell with a radius of 22 pc surrounding WR65 and WR67. The 
infrared Spitzer IRAC image of WR 65 at 8 yum is shown in Fig.[T] 
Despite its unfavorable location close to the edge of the detector, a 
rather compact nebula around WR65 is visible. The radius of the 
nebula is~ 1 .5' or 1 .4 pc. Given the strong mass-loss rate of WR 65 
and its strong UV field, it is plausible that the nebula is physically 
associated with this star. 

In the X-ray sky, WR 65 is located in a very interesting neigh- 
bourhood (see FigO. The star is I'A away from the supernova rem- 
nant (SNR) G320.4-1.2 and 4' away from the X-ray bright "Cir 
Pulsar" and its associated pulsar wind nebula (PWN ). The SNR and 
the C ir pulsar are at the distance d = 5.2 ± 1.4 kpc jGaensler et al.l 
l2002h . This region of the sky is often observed, and serendipitous 
X-ray observations of WR65 are in archives. 

The most homogeneous data set found in archives is from 
Chandra ACIS-I observations conducted from year 2000 to 2005. 
In 2005 WR65 was also observed on one occasion by Chandra 
HRC-I and twice by XMM-Newton, but was not detected during 
XMM-Newton observations affected by high radiative background. 

The local X-ray background is high because WR65 is pro- 
jected on the peripheral region of the PWN mentioned above. 
Therefore we paid special attention to the background subtrac- 
tion .Conveniently, the X-ray emission from the PWN is we l l stud- 
ied jGaensler et"ai]|2002l : I Yatsu et al ] |2005h . lDeLanev et all ( [20061) 
find that it is constant over at least twelve years. We defined back- 
ground regions as annuli around a stellar point source. 

Figure[3] shows the background- subtracted X-ray light-curve 
of WR65 that displays a surprisingly high level of variability. Such 
strong X-ray variability in massive stars is observed only in binary 
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Figure 3. X-ray (0.4-1 0.0 keV) light-curve of WR65. The red dots show 
Chandra ACIS-I count-rates. The blue dots show count-rates converted to 
ACIS-I rates from Chandra HRC-1 and XMM-Newton MOS (upper limit). 



systems l lCorcoran etai]|2005l : lOskinova, Clarke & Pollockll200l[) . 
The WR 65 X-ray light curve gives strong evidence that the star is 
a colliding wind binary (CWB). Unfortunately, the light-curve is 
too sparse to search for the period. Given the time interval between 
the two recorded maxima, one may conclude that the period is not 
longer than about five years. 



3 WR 65 : X-RAY SPECTRAL VARIABILITY 

Five of the Chandra observations of WR 65 yielded enough counts 
for a crude X-ray spectral analysis. The background- subtracted 
spectra obtained at different epochs are displayed in Fig.|4] The 
spectral energy distribution is ha rder than for single WR stars 
jlgnace. Oskinova & Browrill2003l) . The unresolved emission fea- 
tures coincide with the location of lines of Fe, Ca, Ar, S, Si, and 
Mg (cf. Figs|4l[5]l in collisionally ionized X-ray spectra. Interest- 
ingly, the iron lines at xi 6.7 keV (1.8 A) are weaker at some epochs 
compared to others (e.g. they seem to disappear completely on 2005 
February 07). If real, this may reflect the change of the temperature 
in the hot plasma. 

The X-ray spectrum of WR 65 in high states displays strong 
emission features in the 2 -4 A range (see top panel in Fig.|4]l. 
Weak lines of Ca xix {A 3.2A), Ca xx (A 3.0A), and Ar xviii (A 3.1 A) 
are sometimes obser ved in the X-ray spectra of CWBs, e.g. the 
WC8 binary y^Vel jSchild et al."2004f), the WN6 binary WR25 



Raassen. van der Hucht & Mewe.2003j) . and the LBV binary rjCai 
Hamaguchi et al I l2007h . seems plausible to assume that the un- 
resolved Ca and Ar lines contribute to the emission "bump" at 2 - 
4 A in the WR65 spectrum. However, the strength of this emis- 
sion complex in WR 65 is much higher than observed in any other 
CWB. We rule out its origin from the background PWN emission, 
because there are no emission features at 2 - 4 A in the spectra ex- 
tracted from an annulus region around WR65. 

In colliding wind binaries consisting of an OB star and a WC 
star, the latter has a significant ly denser, C, O-enric hed wind that 
is basically opaque for X-rays jOskinova et al.ll2003h . When in the 
course of the orbital motion the wind collision region (WCR) is at 
superior conjunction relative to the WC star (with respect to the ob- 
server), the X-rays generated in the WCR suffer strong absorption. 
Despite many still unclear details, this general picture is observa- 
tionally confirmed by the X-ray spectral analysis of the WC+O type 
binaries WR 140 and Vel. In these systems the strongest photo- 
attenuation is seen in low X-ray states, while the lowest photo- 
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Figure 4. Chandra ACIS-I background subtracted spectra of WR 65 at dif- 
ferent epochs. The date of the observation and the exposure time are given 
in each panel. 



attenuation is found when the X-ray fl ux is high ( iRauw et al.l200(]| : 
ISchild et al.|[2004l : |Pollock et al.ll2005h . 

To investigate how X-ray spectral variations are coupled with 
the variation of the absorption column, we fit a model to the ob- 
served spectra and check how the model parameters vary between 
the different epochs. 

Following a method used in IPollock et all j2005h to suc- 
cessfully analize X-ray spectra of WR 140, we tried to fit the 
observed WR65 spectra alternatively with non-equilibrium ion- 
ization, pshock, and power-law models using the xspec software 
( lArnaud|[T99Q) . However, no meaningful constraints on model pa- 
rameters could be obta i ned b ecause of the poor quality of the 
WR 65 spectra. IZhekovl ( l2007h found that only models with non- 
equilibrium ionization can reproduce the detailed X-ray spectra of 
CWBs. However, he noticed that these spectra can be crudely de- 
scribed by thermal plasma emission with some "average" temper- 
ature. Finally, we select a two-temperature thermal plasma model 
dSmith et alj200lL APEc) where each component is allowed to suffer 
different degrees of attenuation - the "2T2A'h" model. This model 
provides a first insight into the properties of the system by con- 
straining the model parameters. We refer to the plasma with lower 
temperature kTi as the "soft" model component, and to the plasma 
with higher temperature kT2 as the "hard" model component. As 
demonstrated in Figs.|6l|7]the 2T2A'H-model sufficiently well repro- 
duces the shape of the WR 65 spectra observed at different epochs. 

We assume solar composition. Selecting WC-wind abun- 
dances results in scaling down of the parameter that describes nor- 
malization of the emission measure (EM), but the effect of abun- 
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Figure 5. Chandra ACIS-I spectra of WR 65 in high (red) and low (blue) 
states (obtained on 2000/08/14/ and 2004/12/28 respectively). Wavelengths 
of strong spectral lines are indicated. 



dances on spectral shape cannot be discriminated given the poor 
quality of available data. 

The spectrum of WR 65 with highest signal-to-noise and its 
best-fit model are shown in Fig.[6] The temperatures inferred from 
spectral fitting are kTi = 1.0 ± 0.2 keV and kTj = 13.4 ± 7.3 keV. 
Such temperatures can be expected from the collision of metal en- 
riched winds in wide binaries. The temperature is highest at the 

[ oint of head-on collision, and is d ecreasing along the shock cone 
Stevens. Blondin & PoH"oc3ll992h . 

The interstellar absorption column density towards the Cir 
pulsar, A'u.ism = 1 x lO'^^ cm"^, is well constrained and is consis- 
tent with the color excess inferred from our fitting of the WR65 
spectral energy distribution. The neutral hydrogen column densi- 
ties are Nh,i = (2.3 ± 0.3) x 10^' cm"^ for the soft component and 
/V//,2 = (12.2+4.1)xl0^^ cm"^ for the hard component. The column 
density to the soft component, A'h.i, exceeds the interstellar column 
by a factor of more than two, but it is lower than the column density 
to the hard component, A'h,2 . This indicates that the hotter plasma 
is more deeply embedded. 

From the spectrum obtained at 2005/10/18 the model flux is 
Fx(0.6 - 9.0keV) = 1.2 x 10"" erg s"' cm"^, which corresponds to 
Lx = 1.1 X 10^^ erg s"', or in terms of stellar bolometric luminosity 
to logLx/iboi = -7.2 ... - 7.8. The flux of the soft spectral compo- 
nent, Fx(0.6 - 2.0keV) = 8.3 x 10"'^ ergs"' cm"^, corresponds to 
Lx(0.6 - 2.0keV) = 1.0 x 10^' erg s~' which is s imilar to the X-ray 
luminosity of a single O dwarf ( IOskinovall2005l) . 

Spectra observed at other epochs were also fitted with the 
2T2A^H-niodel. For illustration, the spectrum obtained in the low 
state and its corresponding best-fit model are shown in Fig.|7] Fig- 
ure[8] shows the model parameters depending on the time of obser- 
vation. Within the error bars, the temperatures kTi and kT2, and 
the absorption column to the source for the soft component, A'h.i 
are similar at all observed epochs, and therefore are not plotted in 
Fig.[8] In contrast, the column density to the hard component, Nn 2, 
changes at different epochs (panel A of Fig.[8]l. Variations in the 
EM of the hard and soft components are shown in panels B and C 
of Fig.[8] The EM of the soft component changes strongly, reflect- 
ing the X-ray light-curve variations, while the changes in the EM 
of the hard component are less pronounced. 
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Figure 6. ACIS-I spectram of WR65 obtained on 2005/10/18 and the best 
fit spectral model - tbabs(NYi,\)*apec(kT\)+tbabs(N]\2)* apec(kT2). The 
XSPEC spectral fitting package was used. For the model pai'ameters see text. 
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Figure 7. Same as in Fig.|6]for the spectrum obtained on 2004/12/28. 
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Figure 8. Upper panel same as in Fig.O Panel A: variations of absolution 
column density of the hard component obtained from spectral fits. Panel 
B: normalization factor that characterizes the emission measure of the hard 
component. Panel C: normalization factor of the soft component. 



4 DISCUSSION 



The main findings from the analysis of WR 65 X-ray spectra ob- 
tained at different epochs are: (i) the spectra display emission fea- 
tures at the wavelengths of the lines of strongly ionized Fe, Ca, Ar, 
S, Si, and Mg; (ii) the star shows strong spectral X-ray variability; 
(iii) the spectra can be fitted with a two-temperature model, consist- 
ing of a soft component with 0.3;fcri;1.5 keV and a hard compo- 
nent with 2;fcr2;20 keV; (iv) the absorption for the hard component 
is strong and variable, while the absorption for the soft component 
is less severe and nearly constant; (v) the changes in the EM are 
stronger for the soft component. 

Can these observa tions be explained in t he framework of a 
colliding wind model? ICorcoran et alj ( l2005h compare the 2.0- 
10.0 keV light curves of WR 140 and j] Car (both are systems with 
high ec centricity and high inclination) and note that they are very 
similar. IPiTtard et all ( [l998h explain the fight-curve (2.0- 10.0 keV) 
of T] Car using a hydrodynamic simulation of colliding winds. They 
show that the data are broadly consistent with the X-ray emission 
being moderately low throughout most of the orbit. However, a 
sharp peak occurs near periastron, when a greater fraction of the 
stellar winds collides, and the higher pre-shock densities and the 
lower pre-shock velocities result in stronger radiative cooling. The 
eclipse occurs immediately after periastron passage when the com- 
panion of T] Car swings behind the LBV star with its dense wind. 

It is hard to see how a similar model can be used to explain 
the X-ray variability of WR65. Strong variations in the absorbing 
column indicate that the orbit of WR 65 has high inclination too. 
Considering the light-curve (Fig.[3]l we can rule out that WR 65 has 
a period of a few days: the two data points obtained in April 2003 



indicate constant count rates over at least a week. Moreover, two 
data points in the light-curve separated by two months (12/2004 
and 02/2005) are similar, while there is an increase of emission 
by nearly 70% in the next two months. Therefore, WR 65 has most 
likely an orbital period of a few years. In this case the long duration 
of the minimum is inconsistent with its interpretation as attenuation 
by the wind of the WC star immediately after periastron passage. 
Quite different from the short-duration rise of column density ob- 
served in WR 140 and 77 Car, the column density in WR65 is per- 
sistently very high, A'n.ailO"' cm"', sharply decreasing only at the 
maximum X-ray light. 

The temporary variations of column density in W R 65 are sim- 
ilar to those measured in the WC8+0 binary y- Vel l lRauw et"al] 
l200(]|). Its hard X-ray emission increases by a factor of ~ 4 
when column density decreases by factor of ~few during the 
passage of we ak-wind inner side of the sho c k cone across the 
line of sight jWillis. Schild & SteveiisI 1 19951 ; iRauw et al.l I2OOOI; 
ISchild et al.ll2004l) . 

In WR65 the soft X-ray emission is present at all observed 
epochs. Its variation can be explained by changing emission mea- 
sure, while the absorption is nearly constant and comparably low. A 
soft spectral component in the low state is also observed in 7* Vel, 
V444 Cyg and Car. The origin of soft X-rays in these stars is at- 
tributed either to individual stellar winds or to a circumstellar neb- 
ula. Accordingly, in these objects the absorbing column to the soft 
component is similar to the interstellar absorbing column. In con- 
trast, we find for WR 65 that A'n.ism < ^h.i < A'h,2, and the EM of 
the soft component changes with the X-ray light-curve. Thus it is 
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unlikely that the intrinsic wind emission of a companion is solely 
responsible for the soft X-rays from WR 65. 

WR 65 i s a pe rsistent dust maker jvan der Huch 3 l200lh . 
iMonnier et al] ( |2002|) notice that binary stars that are persistent 
dust makers have nearly circular orbits. We can assume that the 
same holds for WR 65. There is a growing number of WC9d binary 
stars that constitute the "pinwhee l" class of objects jTu thill et al.l 
l2008l : lMarchenko & Moffall2007h . In these binary systems the dust 
is formed in the wind collision region and is distributed along a 
wide (compared to the orbital separation) Archimedean spiral. The 
peak of dust production occurs at some distance downstream from 
the colliding wind stagnation point. 

So far there is no direct evidence that WR65 is a pin- 
wheel. However, the topology and structure of circumstellar mat- 
ter in a pi nwheel could help to exp l ain X -ray spectral variability 
in WR65. Iwilms. Allen & McCravl ( I2OO0I) presented a model for 
the absorption of X-rays that includes a treatment of dust grains. 
Their work confirmed that large-grain absorption for soft X-rays 
(< 2keV) is reduced due to the self-blanketing of grains. Interest- 
ingly, from the point of view of radiative transfer this is identical 
to th e reduction of opacity for the X-rays in clumped stellar winds 
jpeld meier. Oskinova & Hamann 2OO30. Winds of WC stars are ex- 
tremely rich in C and O and could have enhanced abundance of 
Mg and Ne. C and O are dominant absorbers for the softer X-rays 
(< 1 keV), while Mg and Ne strongly contribute to the absorption 
of harder X-rays (> 1 keV). In the dusty shell around a WC9d star, 
C would be depleted from the gaseous phase and form grains, re- 
ducing the opacity for the softer X-rays. This may help to explain 
the large differences between the column densities that we deter- 
mine from soft and hard X-rays. 

Recently, iLemaster. Stone & Gardine3 ( l2007h presented 3-D 
hydrodynamical simulations of colliding winds. Scaling the results 
of their simulation to the parameters of WR 65, soft X-ray emission 
could originate from a large volume extending above the orbital 
plane. This emission would show little variation with orbital phase 
and suffer comparably mild absorption in the outer regions of the 
WC star wind, similarly to what is observed in WR 65. On the other 
hand, a high postshock temperature of kT x lOkeV is predicted 
for a confined region near the point where the winds collide head- 
on. Densities up to pJlO"'"* gcm"^^ can be expected in the dusty 
spiral arm s. Assuming the wind opacity for hard X-rays as /<• a; 
5 cm' g ' jAntokhin. Owocki & Brownll2004) . a spiral arm width 
of 10'"* cm would explain the observed high absorption. This width 
is consis tent with the obser ved width of the dusty spiral arms in 
WR 104 l lTuthill et aklboOSl) . The lowest absorbing column can be 
expected at the phase when the line of sight crosses only the outer 
low density part of the spiral arms. It is possible that the Chandra 
observations of WR 65 at high state have caught this orbital phase. 

To conclude, our first brief study of the variable X-ray emis- 
sion from WR 65 indicates that it is a colliding wind system, where 
the dust and complex geometry of its colliding wind region are 
pivotal in explaining its X-ray properties. At present, WR65 is 
the only known WC9d star where X-ray emission is detected. Fu- 
ture observations are needed to constrain the orbital parameters of 
WR 65 and allow for detailed modeling. 
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